SNOTEL Seasonality

Whiskey creek full record

## Warning: Removed 18659 rows containing missing values (geom_path).

Look at raw data

## Warning: Removed 35 rows containing missing values (geom_path).

## Warning: Removed 16 rows containing missing values (geom_path).

How much data is missing?

** A lot of data missing in April let’s throw that out

## Warning: Removed 65 rows containing non-finite values (stat_boxplot).

Look at seasonal monthly median

## Warning: Factor `month` contains implicit NA, consider using
## `forcats::fct_explicit_na`

Seasonal monthly mean

## Warning: Factor `month` contains implicit NA, consider using
## `forcats::fct_explicit_na`

look at max each water year

## Warning: Removed 3 rows containing missing values (position_stack).

Beaver Spring and WHiskey Creek seasonality

Beaver Spring

PRISM, snow course, snodas, snotel

## Warning: position_dodge requires non-overlapping x intervals
## Warning: Removed 18027 rows containing missing values (geom_col).
## Warning: Removed 2192 rows containing missing values (geom_path).

plot PRISM with Snow course SWE

## Warning: Removed 14206 rows containing missing values (position_stack).

## Scale for 'y' is already present. Adding another scale for 'y', which
## will replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which
## will replace the existing scale.
## Warning: Removed 3 rows containing missing values (geom_point).

plot SNOTEL SWE and PRISM and SNODAS

# SNOTEL and PRSIM
bs_snotel_prism <- bs_sc_snotel_prism_snodas %>% 
  mutate(month = month(date, label = TRUE, abbr = TRUE)) %>% 
  filter(month %in% c("Nov", "Dec", "Jan", "Feb", "Mar")) %>% 
  select(date, month, SNOTEL, PRISM) %>% 
  drop_na() %>% 
  mutate(Date = date) %>% 
  addWaterYear()

graph_correlation(bs_snotel_prism, x_data = bs_snotel_prism$PRISM, 
                        y_data = bs_snotel_prism$SNOTEL,
                        title = "Beaver Spring SNOTEL Daily SWE and PRISM ", 
                         xlabel = "PRISM Accumulated Snowfall (mm)",
                         ylabel = "SNOTEL SWE (mm)") +
  scale_y_continuous(expand = c(0,0), limits = c(0,400))+
  scale_x_continuous(expand = c(0,0), limits = c(0,400))+
  facet_wrap(~waterYear)
## Scale for 'y' is already present. Adding another scale for 'y', which
## will replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which
## will replace the existing scale.
## Warning: Removed 45 rows containing missing values (geom_point).

# SNOTEL and SNODAS
bs_snotel_snodas <- bs_sc_snotel_prism_snodas %>% 
  mutate(month = month(date, label = TRUE, abbr = TRUE)) %>% 
  filter(month %in% c("Nov", "Dec", "Jan", "Feb", "Mar")) %>% 
  select(date, month, SNOTEL, SNODAS) %>% 
  drop_na() %>% 
  mutate(Date = date) %>% 
  addWaterYear()

graph_correlation(bs_snotel_snodas, x_data = bs_snotel_snodas$SNODAS, 
                        y_data = bs_snotel_snodas$SNOTEL,
                        title = "Beaver Spring SNOTEL Daily SWE and SNODAS ", 
                         xlabel = "SNODAS SWE (mm)",
                         ylabel = "SNOTEL SWE (mm)") +
  scale_y_continuous(expand = c(0,0), limits = c(0,500))+
  scale_x_continuous(expand = c(0,0), limits = c(0,500))+
  facet_wrap(~waterYear)
## Scale for 'y' is already present. Adding another scale for 'y', which
## will replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which
## will replace the existing scale.
## Warning: Removed 20 rows containing missing values (geom_point).

compare max SWE snow course to PRISM max accumulated

## Error in as.character(x): cannot coerce type 'closure' to vector of type 'character'

Whiskey Creek

plot PRISM with Snow course SWE

## Warning: Removed 14155 rows containing missing values (position_stack).

## Warning: position_dodge requires non-overlapping x intervals
## Warning: Removed 17976 rows containing missing values (geom_col).
## Warning: Removed 2051 rows containing missing values (geom_path).

## Scale for 'x' is already present. Adding another scale for 'x', which
## will replace the existing scale.
## Scale for 'y' is already present. Adding another scale for 'y', which
## will replace the existing scale.
## Warning: Removed 2 rows containing missing values (geom_point).

## Scale for 'colour' is already present. Adding another scale for
## 'colour', which will replace the existing scale.

compare SWE snotel to swe snow course for whiskey creek and beaver springs

## Scale for 'x' is already present. Adding another scale for 'x', which
## will replace the existing scale.

## Scale for 'x' is already present. Adding another scale for 'x', which
## will replace the existing scale.
## Scale for 'y' is already present. Adding another scale for 'y', which
## will replace the existing scale.

## Scale for 'y' is already present. Adding another scale for 'y', which
## will replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which
## will replace the existing scale.

## Scale for 'y' is already present. Adding another scale for 'y', which
## will replace the existing scale.

compare SWE SNODAS to swe snow course for whiskey creek and beaver springs

## Scale for 'y' is already present. Adding another scale for 'y', which
## will replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which
## will replace the existing scale.

## Scale for 'y' is already present. Adding another scale for 'y', which
## will replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which
## will replace the existing scale.
## Warning: Removed 2 rows containing missing values (geom_point).

Fluted Rock

## Scale for 'y' is already present. Adding another scale for 'y', which
## will replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which
## will replace the existing scale.

## Scale for 'y' is already present. Adding another scale for 'y', which
## will replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which
## will replace the existing scale.

## Scale for 'y' is already present. Adding another scale for 'y', which
## will replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which
## will replace the existing scale.

## Warning: position_dodge requires non-overlapping x intervals
## Warning: Removed 5561 rows containing missing values (geom_col).

Bowl Canyon

## Scale for 'y' is already present. Adding another scale for 'y', which
## will replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which
## will replace the existing scale.
## Warning: Removed 17 rows containing missing values (geom_point).

## Warning: position_dodge requires non-overlapping x intervals
## Warning: Removed 5558 rows containing missing values (geom_col).